Squeezed light from spin squeezed atoms 
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We propose to produce pulses of strongly squeezed light by Raman scattering of a strong 
laser pulse on a spin squeezed atomic sample. We prove that the emission is restricted to a 
single field mode which perfectly inherits the quantum correlations of the atomic system. 
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Squeezed light and entangled beams of light can be 
used to probe matter and to study mechanical motion 
with better resolution than classical light Entan- 
gled photon sources can be used for lithography with a 
resolution below the optical wavelength ||, ^, and the 
active field of quantum information profits from the de- 
velopment of non-classical light sources [§, ^, ^- Non- 
linear crystals in optical parameter oscillators and laser 
diodes with suitable feed-back have been the work horses 
in most experiments on non-classical light. As a figure 
of merit for the degree of non-classicality of these light 
sources, one may refer to the noise suppression observed 
in direct or homodyne photon detection measurements. 
Compared to classical sources it has so far been possi- 
ble to reduce the noise (variance) by about one order of 
magnitude. In order to make a significant difference in 
practical applications, further noise reduction is really 
necessary. 

It has been proposed that the large non-linearity and 
low absorption in resonant Raman systems can lead to 
ideal four- wave mixing and substantial squeezing 0, ||, §| . 
As an alternative approach we propose in this Letter 
to use spin squeezed ensembles of atoms as sources of 
squeezed light. Atoms can be entangled in such a way 
that the fiuctuations in occupancies of different internal 
states are significantly suppressed. This phenomenon is 
refered to as spin squeezing, because a two level atom can 
be described formally as a spin 1 /2 particle, and the inter- 
est in spin squeezed states arose already a long time ago 
in connection with ultra-precise spectroscopy and atomic 
clocks . Squeezing of spins was originally believed to 
be very complicated, but recent proposals based on quan- 
tum non-demolition measurements of atomic populations 
| pj| , on coherent interactions in Bose-Einstein conden- 
sates [ p^ p^ , and on interactions between laser excited 
atoms []14| have changed this impression and suggested 
that really significant spin squeezing is achievable. The 
main purpose of the Letter is to demonstrate that the 
atomic quantum correlations can be perfectly transferred 
to the field. This result is readily obtained within a sim- 
plified model where both atoms and field are described 
by single harmonic oscillators ||l5|] but we show that cor- 
relations in the atoms can be mapped perfectly on the 



field also in an a priori multimode situation. 

The emission of light is treated by a simple gener- 
alization of the theory of stimulated Raman scattering 
Ip^ , [l7| , |l8|| to and this part of our proposal can be 
analyzed without specifying the model for spin squeez- 
ing. Our ensemble of two-state atoms is assumed to be 
strongly elongated and it is treated in a ID approxima- 
tion. It is illuminated by a strong laser field Eg propa- 
gating along the z-axis of the system. This opens up a 
channel for an atom in the b state to go to the a state by 
absorbing a photon of frequency lUs and emitting a pho- 



ton of frequency ojq 



As a result, a field at this 



frequency builds up and propagates through the sample. 
This process is described by the following coupled set of 
equations for the atomic and field operators 
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Eq{z,t) = -iK2Wb{z,t)ll>a{z,t)Es{z,t). 



(2) 



where m = J2il''aif'bi/{fi'^'!^i) and K2 = 2TrhuqKi/c. fij^ 
are dipole moments of the atomic transitions and are 
the (large) detunings with respect to intermediate levels, 
see Fig. 0. i}a{z,t), i)l{z,t) and ^pb{z,t), i'l{z,t) are 
annihilation and creation operators for atoms in states a 
and b; Wb{z^ t)ipa{z, t) is the positive frequency part of the 
atomic dipole operator, taking into account the atomic 
density at position z in the ensemble. In Eq. (|^) we have 
assumed that there is no dephasing of the ab coherence 
and we have assumed the validity of the slowly varying 
envelope approximation for the emitted field in Eq. (^). 

We restrict our analysis to the case where the atoms 
are almost entirely in the a state when the Eg field is 
applied. This implies that the population difference ap- 
pearing in Eq. (|l|) can be replaced by the density of 
atoms n(z). This density we represent as a c- number 
throughout the duration of the output coupling. This 
allows us to define a dipole operator "per atom" by 
ipl{z,t)ipa{z,t) — n{z)Q{z,t) and we obtain the linear 
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operator equations 
d 



dt 



Q{z,t) 



d_ Id 

dz c t 



-zKlE:{z,t)E,{z,t) (3) 
-iK2n{z)E,{z,t)Q{z,t). (4) 



If we generalize the analysis in |lq, 17; 18 to inhomo- 
geneous media, we can solve Eqs. (S^^and ^ analytically 
in terms of the input field Eq{0, t) at the entrance face of 
the sample at z = and the initial position dependent 
atomic polarization Q{z,0). In particular, for the field 
operator we get 



Egiz,T) = Eg{0,T)-lK2E,{T) I Vl(^',0)V'a(^',0) X 



X Jo \ 2^a{T) J n{z")dz"j dz' (5) 

where the new time coordinate is r = t — z/c, Jo(') 
is a Bessel function of the first kind, and a(T) = 

KlK2j^\E,iT')\'dr'. 

The expression (|^) must be evaluated at the position 
z — L of SL detector outside the atomic sample. In the 
absence of atoms, the field equals the incident quantum 
vacuum field i?g(0, r). The atomic sample is able to re- 
place the vacuum with an entirely different field. In order 
to analyze the quantum properties of Eg it is convenient 
to imagine a time integrated homodyne dection at the 
detector. By choosing the temporal form of the strong 
local oscillator field in this detection we select a certain 
spatio-temporal mode of the field represented by the field 
operator 



£*{T)EgiL,T)dT 



(6) 



where / |f pdr — 1. In choosing £{t) we should seek 
to ensure that a is a mapping of the precise collec- 
tive operator of the atomic sample that is known to be 
squeezed. Mathematically, it is easy to show that in or- 
der to probe Jh{z')Q{z' ,0)dz' we should choose £{t) as 
the normalized solution of Eqs. with the initial con- 
dition n{z')Q{z', 0) replaced by h{z'). Physically, this re- 
flects that the question of mode matching coincides with 
the problem of identifying the classical field radiated by 
a classical dipole distribution. In particular, if the uni- 
form integral of atomic operators / ipl{z' ,0)ipa{z' ,0)dz' 
is squeezed, mode matching is accomplished by taking 
£{t) to be the (normalized) solution to the classical Ra- 
man scattering problem. In that case we get the mapping 



J- 



(7) 



where J- = J 'il'l{z, t)il'a{z, t)dz, and N is the total num- 
ber of atoms. 



It is a natural concern, whether a 3D analysis will pre- 
serve the possibility to couple to only a single spatio- 
temporal field mode. The issue has been addressed in the 
classical case, where a diagonalization of the first order 
coherence function (Karhunen-Loeve transformation) of 
the field indeed shows that a single mode dominates the 
output, provided the Fresnel number of the (active part of 
the) atomic sample is of the order of unity Q . Choosing 
the experimental parameters accordingly we thus expect 
this result to apply also for the quantum field in three 
dimensions. 

We recall that several schemes now exist for the gener- 
ation of non-classical states of atomic spins. The above 
analytical treatment is general and independent of the 
method of spin squeezing, and as shown by Eq.(7) (valid 
only if the majority of atoms occupy the state a), the 
atomic state is transfered perfectly to the light field. 
Mean values and variances for the field observables are 
therefore explicitly known, and, e.g., the orders of mag- 
nitude squeezing derived in |l^ apply to the emitted 
light pulse. It may be useful, however, to relate the emit- 
ted quantum field directly to the dynamical variables in 
the spin squeezing process. This may provide further in- 
sight in the origin of squeezing; it may provide a useful 
tool for the application of the squeezed light as input to 
another quantum system |19, and, it may be useful 



to analyze processes where the spin squeezing and the 
light emission occurs simultaneously. 

As an example, consider spin squeezing by coUisional 
interactions in a two-component Bose Einstein conden- 
sate First, a Bose Einstein condensate is formed 
in only one of the internal states a. By a short reso- 
nant Raman pulse, all atoms are transfered to an equal 
superposition of a and b. We assume that the inter- 
action strengths between the atoms are not all equal, 
daa = 9bb 7^ .9ah, SO that the interaction terms in the fully 
quantized interaction Hamiltonian 



H 



■ 9ab4'l{r)4ir)Mr)Mr)} (8) 



cause a Kerr-like phase evolution of amplitudes on states 
with different numbers of atoms in the two internal 
states. This results in squeezing of an appropriate col- 
lective spin variable as pointed out by S0rensen et al. 

0- . . . 

As in our earlier work on this problem |13| we use the 

positive P method [ pl| |2^ , p3| to describe the squeez- 
ing process. This means that averages of all normal 
ordered operator products can be calculated as ensem- 
ble averages by the substitution of atomic field operators 
xjj = (■(/'a, 'i/'J, ■;/){,, 'i/jj) by pairs of two-component "wave- 
functions" , xj} = {ipai , V'a2 , V'fci 7 "062 ) . The dynamics of the 
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wave functions is given by four coupled and noisy "Gross- 
Pitaevski" equations, see details in Ref. Starting 
from a coherent initial state of all atoms in an equal su- 
perposition of internal states a and b, we numerically 
simulate solutions of these equations to obtain an en- 
semble of i/j's describing exactly (up to sampling errors) 
the quantum correlations of the system. 

Once a sizable spin squeezing is obtained we want to 
transfer this special quantum state to a light pulse, and 
to be able to use the results of the above analysis we 
first bring the internal state of the atoms close to the a 
state, i.e., with a new resonant Raman pulse we rotate 
the collective spin close to the north pole of the Bloch 
sphere. This rotation is applied to the individual sets of 
"wavef unction" realisations of the simulation {i — 1,2): 

with e = 7r/2. 

The atom field interaction in the coupled equa- 
tions (|^) and (||) leads to a natural positive P represen- 
tation of the field Eq{z,t): Replace in Eq. ^ Eq{z,t) 
by a c-number field Eqi{z,t) and -0^(2;', O)'0a(2:', 0) by 
'>Pb2iz' ,0)'ipai{z' ,0), and make a similar replacement in 
the hermitian conjugate equation where i/Jaiz' , 0)ipb{z' , 0) 
is replaced by ipa2{z' ,0)'4'bi{z' ,0) to yield the c-number 
field Eq2{z,t). In both cases, the incident vacuum field 
can be represented by a zero, since the positive P repre- 
sentation yields normally ordered expectation values as 
simple products. From our original ensemble of quadru- 
ples of wavefunctions xp describing the atoms just before 
we turn on E^ we obtain in this way an ensemble of pairs 
{Eqi{z,T), Eq2{z,T)) describing the generated light field 
Eq{z, t). With this ensemble any normally ordered field 
expectation value can be found. 

As an example we have investigated 2000 atoms of 
mass m and with ID interaction strengths {gaa, gab, gbt = 
(1.0,0.5, 1.0) X 5 X lO^'^hfllo where fl is the frequency of 
the harmonic trap and Iq — y^h/mfl is the associated 
characteristic length. The atoms are spin squeezed by 
collisional interactions for a time t = 3.051^^. They are 
subsequently driven towards the state a, and hereafter 
they are illuminated with the Eg light which builds up 
a maximum strength of E^ax — y^2TThuj / lo in a time 
of roughly trise — IOO/q/c. The matter- light coupling is 
chosen to be ki = 10~'^c/27r/ia;. 

As we know the atoms to be spin squeezed we also 
know one mode of the Eq field which will be squeezed: 
The one corresponding to the simple uniform integral 
of J ipl{z' ,0)il)aiz' ,0)dz' as described in the discussion 
after Eq. (||). Even with the simplifications leading to 
Eqs. (H) and (||) the atomic state could in principle 
radiate into many other modes and do so with vary- 
ing quantum statistics. To check whether such other 
modes are present in this case we calculate the first or- 
der correlation function of the field {E'^{L, T')Eq{L, r)) = 



Eq2{L,T2)Eqi{L,Ti) ( . . . iudicatcs averaging over the 
positive P ensemble). When it is diagonalized we find 
that almost all population is in fact in the expected mode 
£(r) which is plotted in Fig. [l[ In other scenarios it might 
be difficult to calculate beforehand exactly which collec- 
tive atomic operator is squeezed and then an analysis 
like this would be necessary in order to pick the local 
oscillator field for the homodyne detection. 

Having confirmed that only one mode is populated we 
now turn to the quantum character of the field. It de- 
pends on how we choose in Eq. (|): For 6* = 7r/2 it 
will approximate squeezed vacuum, for slightly differ- 
ent from 7r/2 it will approximate a squeezed coherent 
state. As described in [|l^ the spin squeezing ellipse is at 
an angle to the coordinate axes. This carries over to the 
light-field which is squeezed in an appropriate quadra- 
ture component = (ae*"* + a^e'^'^)! ^/2. With the 
above parameters, we find from the positive P simula- 
tions that the minimum variance is 0.04 ± 0.01 corre- 
sponding to a reduction by a factor of more than 12 from 
the standard quantum limit. To illustrate the positive P 
results in Fig. ^ we show both a histogram and a scat- 
ter plot obtained from the positive P ensemble of pairs 
(ai,a2) representing (a, a^^) for the mode of the field. a\ 
and ail on average, are complex conjugate quantities, so 
that the expectation values of hermitian field operators 
(AT^) — aie**^ -I- 026"*"^ is real for all 0. To represent 
our simulated results we have made a histogram for the 
values of a; = Aq and p — X^ji, indicated by the real 
part of (fli -+ / \/2 and {iai — 102) /%/2. The histogram 
in Fig. 2 has been obtained from 10^ independent re- 
alizations of the noisy Gross-Pitaevskii equation for the 
atomic spin squeezing, while the scatter plot contains 
only 6800 points each representing a single realization. 

The perfect output coupling of atomic correlations is 
the main result of the paper, and we imagine that many 
other atomic quantum states may be taken as start- 
ing point for non-classical light generation, e.g., EPR- 
correlated separated atomic ensembles, and that the 
mechanism may be used also for reliable interspecies 
teleportation. A squeezed light pulse may also be ap- 
plied at the "dark" input port of a beam splitter, to 
cause the splitting of a strong pulse (with the appropriate 
mode function £{t)), into two twin pulses with perfectly 
matched photon statistics. It was recently demonstrated 
that a classical light pulse can be brought to a complete 
halt in an atomic sample and that it can be subsequently 
released |Q. Our study suggests that a quantum light 
pulse can be similarly stored and retrieved, and that one 
may perhaps introduce a known or unknown signal field 
in a condensed sample, manipulate its quantum state 
by collisional interactions and release it with the orig- 
inal mode function preserved. We emphasize that the 
spin squeezing mechanism is not crucial for the squeezed 
light output, and for example the QND-atomic detection 
and spin squeezing scheme of Kuzmich ct al ||Tl|] may be 
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an interesting possibility. Some schemes for atomic spin 
squeezing may even be compatible with a transmission 
or production of light with a group velocity so low, that 
it may constitute a cw source of squeezed light. 
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FIG. 1: The shape of the Es pulse and the mode function 
£ of the emitted Eq pulse. The mode function is found by 
diagonahzing the first-order correlation function of the field, 
but coincides with the signal expected when the single atom 
operator Q is uniform over the sample. The insert shows the 
atomic level scheme of the proposal. 




FIG. 2: Histogram illustrating the prediction for the squeez- 
ing of light. The quadrature components of the output 
field axe represented by pairs of numbers Re(ai -|- a2)/%/2, 
Re(mi — ia-i)/\/2, the distribution of which forms a squeezed 
ellipsoid shape in phase space. The actual amount of squeez- 
ing cannot be directly determined from this plot; it requires a 
computation of mean values and variances, making use of the 
fact that the distribution of the complex numbers a\ and a2 
represent mean values of normally ordered field operators. 10^ 
realizations contributed to the histogram and in the insert is 
for illustration shown a scatter plot of 6800 of the representive 
points. 



